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Scattering is an important phenomenon which is observed in systems ranging from the micro to 

Lfi, ] macro scale. In the context of nuclear reaction theory the Heidelberg approach was proposed and 

later demonstrated to be applicable to many chaotic scattering systems. To model the universal 

properties, stochasticity is introduced to the scattering matrix on the level of the Hamiltonian 

^sj ' by using random matrices. A long-standing problem was the computation of the distribution of 

the off-diagonal scattering-matrix elements. We report here an exact solution to this problem and 

present analytical results for systems with preserved and with violated time-reversal invariance. Our 

derivation is based on a new variant of the supersymmetry method. We also validate our results 

^^ , with scattering data obtained from experiments with microwave billiards. 
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A large part of our knowledge about quantum systems comes from scattering experiments. Furthermore, even in 
classical wave systems, observables can often be traced back to a scattering process [l|. Important examples stem from 
nuclear, atomic and molecular physics, mesoscopic ballistic devices, microwave and elastomechanical billiards [lHj]. 
Accordingly, the investigation of scattering phenomena has been a subject of major interest from both theoretical and 
experimental points of view [l|-|24j . 

The scattering matrix (S-matrix) is the key quantity in scattering theory. It relates the asymptotic initial and 
final Hilbert spaces spanned by a quantum system undergoing the scattering process. Owing to the flux conservation 
requirement the S'-matrix is unitary, i.e., SS^ — S^S = 1. The universal features of chaotic scattering systems call 
for a statistical description. The Heidelberg approach provides a convenient basis for this. The associated S'-matrix 
elements are obtained in terms of the Hamiltonian H describing the scattering center [5J , and the coupling vectors 
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W c which incorporate the information concerning the interaction between the internal states of H and the M open 
channel states labeled c = 1, ...,M. The full system resides in the latter asymptotically before or after the scattering 
event. The S-matrix elements are obtained as [a, LL2| 



S ab {E) = S ab - i2nW^G(E)W b , (1) 

where the inverse of the propagator G(E) reads 



M 



G- l {E) = E\ N -H + m^W c Wl. (2) 

Tt : c - x 

For the sake of simplicity we restrict the considerations to the case where the average S-matrix is diagonal 25] . This 
scenario translates to the condition that the coupling vectors W c are orthogonal, viz., W\Wd = {lc/^)^cd]c,d = 
1,...,M@,[H. 

Stochasticity is introduced on the level of the Hamiltonian describing the internal dynamics [5]. In view of the 
universality conjecture for quantum chaotic systems [26j . the Hamiltonian H can be modeled by the Gaussian ensemble 
of N x N random matrices with the distribution of H given as [J, [271 [28| , 



V(H)<xex V (~^trH 2 ). (3) 

Here, v 2 fixes the energy scale and the index /3 signifies the appropriate symmetry class, i.e., the invariance properties 
of the Hamiltonian. We focus on the cases [3 = 1 (orthogonal ensemble) and ft = 2 (unitary ensemble) which apply 
to "spinless" systems with, respectively, preserved and violated time-reversal (T) invariance [4|, [271 128( . 

In their pioneering work 6], using the supersymmetry method [29], Verbaarschot et al. calculated the energy 
correlation function of two S-matrix elements. Further progress in this direction was made in [7j where three- and 
four-point S-matrix correlation functions were evaluated. While these provide rich information about the scattering 
process, a full statistical description requires the knowledge of the distributions of the S-matrix elements. The 
complexity involved in the calculations of such energy correlation functions [a, [71, |3y] indicates that the derivation of 



the S-matrix distribution functions is an even more challenging task. However, it was partially accomplished in 13 1 



where distribution of the diagonal S-matrix elements was derived. Moreover, in 14| the statistics of transmitted 
power which is given by \G nm (E)\ 2 ,n ^ to, was calculated. These results as well as those for the correlation functions 
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in systems with preserved and violated 7" invariance have been verified in microwave experiments 

unsolved problem was the computation of the distribution of the off-diagonal S'-matrix elements which could not 

be tackled with the well-established methods of Refs. [y, |7j, [l3j. In the present work we provide analytical results 

for those of their real and imaginary parts. The novelty of our approach lies in that a nonlinear sigma model is 

constructed based on the characteristic function associated with the distributions which is the generating function for 

the moments. In contrast, the standard supersymmetry approach starts from the generating function for the S'-matrix 

correlations. 

We introduce the notation p s (S a b), with s — 1, 2 to refer to the real and imaginary parts of S a b, respectively. Thus 
Eq. ([1]) yields for the off-diagonal (a ^ b) elements 



Ps(Sab) = n((-i) 3 W*GW b + i s Wl&W a ). (4) 

which we denote by P s (x s ), involves the nontrivial task of performing an 



Determining distributions for p s (S a b 
ensemble average, 



P s (x s )= J d[H}V(H)S(x s -p s (S ab )), 8 = 1,2. 
We instead first compute the corresponding characteristic function, 



R s {k)= l d[H]V(H)ex V (-ikp s (S ab )), 

and then obtain P s (x s ) as the Fourier transform of R s (k), 

1 f°° 
P s (x s ) = — / dkR s (k)exp(ikx s ). 
2tt J_ 00 

Defining the 27V-component vector W and the 2N x 2N matrix A s as 



W 



W a 
W b 



A„ = 



{-i) s G 
i s G^ 



we rewrite R s (k) as 



R s (k)= I ' d[H]P(H)&q?(-ikwW i A a W). 



(5) 



(6) 



(7) 



(8) 



(9) 



In order to perform the ensemble average, we map the statistical model to superspace. For this we introduce the 2JV- 
vectors z T = \z^,z^\ and £ T = [CjiCif] consisting of complex commuting and anticommuting (Grassmann) variables, 
respectively. The supervector is constructed in the usual manner 3l| as ^ T = [z T , £ T ]. With the aid of these vectors, 
and using the standard multivariate Gaussian integral results for commuting and anticommuting variables, we recast 
the characteristic function as 

R s (k) = /d[*]exp(|(W t \E' + * t W)) f d[H]V{H)exp f-L-tft A -i#Y (i ) 

Here Aj 1 = I2 <8> ^J 1 , and W^ = [W',0] is a 4/V-component vector. Generally, the ensemble average can be done 
with the help of the supersymmetry method [y, Il2l |31 |. thereby leading to an enormous reduction in the degrees 
of freedom. However, in the present case the averaging is not straightforward because the matrix A J 1 is not block 
diagonal. To resolve this problem we introduce the transformations 



z^T+z, *t_>*t ; C _,|T-C,C t ^C t , 



T± = 



±(-i) s l N 
~i s t N 



(11) 



Here, we have used the fact that the complex quantities z (£) and z^ ((^) are independent of each other. The Jacobian 
of the transformation is (— \) N 2~ 2N for /3 = 1 and equals (— 1)^ for (3 = 2. After application of Eq. (jlip we have to 
distinguish between the orthogonal and the unitary cases. 



For /3 = 2 we obtain 

R s (k) = (-1) N /d[*]exp(^(Ui* + * t W)) J d[H]V{H)e* v (^-^A- 1 *), (12) 



-*»w 6 t ,H) a w2,o,o 



and matrix 



with the 47V-dimensional vector U| = 

diag [-(G- 1 )t,G- 1 ,-(G- 1 )t,-G- 1 ], and * T = [z T ,( T ] as above. 

For j3 = 1 we decompose the 2iV-vector z into its real and imaginary parts x and y (not to be confused with X\ 
and X2) to construct a 47V-vector. In addition, we symmetrize the vector £ using C*,Cb along with CaiO» thereby 
also doubling its size. The associated Jacobian equals 2 2N and thus cancels that of the transformation Eq. (fTTj) . The 
8iV-dimensional supervector is given as \F 1 ' = [x^, y„ , xj, y^ , £|, — £j, Q , — C^] and 

i? s (fc) = (-l) w fd[V]exp(i&V s ) / , d[H]P( J ff)exp(-^-* t ^- 1 *y (13) 

In this case V^ = [i s wj, -i s+1 w+, w+ , iw~,0, 0,0, 0], where u?± = {{-i) s Wj ± W^ T )/2, and A" 1 = 
diag[-(G" 1 ) t ,G" 1 ,-(G" 1 ) t ,-G" 1 ] ® 1 2 . Since the matrix ^l^ 1 is block diagonal in both Eqs. ([12]), CCU the 
ensemble averaging now is straightforward. 

Next we employ the standard supersymmetry techniques [3l| using the Hubbard-Stratonovitch identity to map 
the integral over the 8iV//3-dimensional vector \P to a matrix integral in superspace involving an 8//3-dimensional 
supermatrix a of appropriate symmetry. This yields 

R s (k)= / d[cr] exp ( - r str cr 2 - - str InS- -f\, (14) 

i M E 

Y, = <j E ®\ N + —L®^rw c Wl, a E = a - —. -l 8 /s, 
4/e *— : ' 4nk ' 



with str denoting the supertrace. Here r = {Aj3-K 2 k 2 N)/v 2 , and L = diag(l, — 1, 1, — 1) eg) 12//3- F s equals 

VjL- 1 /2S-iL- 1 / 2 V s for /3 = 1, and UjL- 1 /2 S -i L -i/2 W for ^ = 2 with L = L ® Ijv- The supersymmetric 
representation, Eq. (fl4| . forms one of the key results in this letter. 
Using the orthogonality of W c we eventually obtain 

M 

str InS = iVstr lno- B + ^2 str ln ( 1 8/ / 3 + t\ ,j ~e Ij 

c=l 

M M 

5T 1 = (jj, 1 ® Hat - a E x ® V —W C W\ + V p (c) <g) — WcWJ, (15) 

with p^ c ' = (<te + i7 c /(47rfc)L) _1 . Furthermore, F s equals a linear combination of matrix elements of p^ multiplied 
with 7c, where c — a, b. Note that in Eq. (fT5")) the first term is of order N while the remaining terms are of order 
1. Thus, in order to perform the limit N — > 00, we may apply the saddle point approximation. This leads to a 
separation of a into Goldstone modes uq and massive modes [3l|. The integrals over the latter, being Gaussian ones, 
can be readily done and yields unity. We are therefore left with an expression involving only the Goldstone modes, 
and consequently our sigma model reads 

„ M 

R s {k) = / d^(a G )e- i F » J] sdet- - (l s/p + TT^L) > (16) 

•* c=l 7F 

with sdet denoting the superdeterminant and a replaced by og in all the ingredients of Eq. (|14p . In order to perform 
the remaining integrals we proceed as in [y, |31| and express oq in terms of an 8//3-dimensional supermatrix Q as 
a G = (E/8Trk)t 8/ p~(A/8Trk)Q with Q 2 = -l 8 /p. Here, A = (4v 2 -E 2 W? with A/(2ttv 2 ) identified as the celebrated 
Wigner semicircle. We use the same parametrization of Q as in [y, [l2J, |32J . For /? = 2, it involves pseudo eigenvalues 
Ai G (l,oo),A2 € (—1,1), angles (f>i,(f>2 S (0,27r) and four Grassmann variables. For /? = 1 we have three pseudo 
eigenvalues Ao € (— 1, l),Ai,A2 € (l,oo), two 0(2) angles </>i,</>2 G (0,27r), three SU(2) variables m,r,s e (—00,00), 



and eight Grassmann variables. The product over the superdeterminants in Eq. (|T6f involves the pseudo eigenvalues 
only, and turns out to equal 

M + + A 

^S^+aoW+a^ for/3=1 ' 

M + I \ 

Here gf = (v 2 ± 7 2 )/(7 c A). g+ is related to the transmission coefficient T c = 1 — 1 5" cc | 2 as g+ = 2/T c — 1 6]. The 
exponential part in Eq. (fT6|) also involves other variables and, in fact, is quite complicated for (3 = 1. 

For (3 — 2 the integrals over the Grassmann variables and over the angles can be performed and we obtain the same 
distribution for the real and imaginary parts of the S'-matrix elements, 

f°° C 1 k 2 / I \ 

R s {k) = l- J dX x J dX 2 _ Jv(\ u A 2 ) (titl + tit 2 ) J (kyjtltl) , (17) 



where J n {z) represents the nth order Bessel function of the first kind, and t J c = \/\Xj — l|/(.9c + Aj), j = 1,2. The 

"1" in Eq. (fTT|) is an Efetov-Wegner contribution [29( which is essential for the correct normalization, R s (0) — 1. The 
distribution function is obtained using Eq. ([7]). It is given as 

P s ( Xs ) = d 2 f( Xs )/dxl, 

f(x)=x@(x) + ["dXi f dX 2 Zn Xu \ 2 l (^ + &S) («& - -T 1/2 ®titl ~ * 2 )- (18) 



-l 



4tt(Ai - A 2 ) 2 



Here, 0(u) is the Heaviside function. Since the distributions of the real and the imaginary parts of the S'-matrix 
elements are identical in this case, it is evident that their phases have a uniform distribution and that the joint density 
of the real and the imaginary parts is a function of \fx\ + x\ only. The last observation can be used to calculate the 
distribution of the modulus [33] . 

For (3 — 1, the calculation involved is rather cumbersome. Nevertheless, we managed to perform all but four 
integrals. For the real part the result reads 

-1 />1 />00 />00 /> 2 7T 4 

i? i( fc ) = 1 + ^/ dX o d\ x d\ 2 d^(A ,A 1 ,A 2 )J'o(Ao,A 1 ,A 2 )J] K „r, (19) 

where 

J= (1- Ag)]Ax- A 2 | 



2(A? - 1)V2(A| - l)i/2(A x - A )2(A 2 - A ) 2 ' 
The k's entering Eq. (fT9|) are functions of 



P 3 c= J +— r , J =0,1,2; p±=^±^, 
8(ffc + Aj) 

+ ifE ■ -\r i i 



VA J " S M 5 + + Ai 5 + + A 2 5+ + A 

and of the complex conjugate of q^, rf = (q^)* , and the quantities u> = 2^/XY, I = X/Y, m = Y/X, where 

X = 2pt + q-e-^+r-e^, 
Y = 2p+-q-e^-r-e-^. 

It can be verified that uj 2 is real and takes values from the interval [0, 1] for all values of the parameters involved. The 
k's are given as 

Ki = «ii Ji(kuj), k 2 = K 2 iJ {kuj) + K 22 J 2 (kco), 
k-3 = K 3 iJi(kuj) + n 32 J 3 (kuj), 

K 4 = K41 Jo(kuj) + K4 2 J 2 (k(jj) + Ki 3 Ji(kljj) 1 



with the entries Kij, i,j = 1,-- • ,4 given in the appendix. 

The characteristic function for the imaginary part, Pt2(k), is obtained by multiplying —i to the right hand side of 
the expressions for qf in Eq. (f2TJ)) , and changing rf accordingly. The distributions of the real and the imaginary parts 
of the S-matrix elements can be cast into the form, 



P s {x s ) = 5{x s ) + - — 



dh , d 2 h , d 3 f 3 , 5 4 / 4 



9a;? <9x? 



ax! 



(21) 



/i = {khx s /uj) , 

/2 = -<«2i + « 22 (l-2^/a; 2 )>, 

/3 = ~([k31 + ^32(3 - 4:X 2 Juj 2 )]x s /uj) , 

,A = ( [ku + k 42 (1 - 2^/w 2 ) + k 43 (1 - 8a; 2 /w 2 + 8xj/uj 4 )]). 



Here the angular brackets represent the following: 

"1 />OC />00 /•27T 



w = 



1 



16tt 



d\ / dAi / rfA 2 / ^J(A ,Ai,A 2 )J- (Ao,A 1 ,A 2 )2^ 2 -a; 2 )- 1 / 2 e(w 2 -a;2). 
1 Ji 7i Jo 



Details of the supersymmetry calculations and further results will be given elsewhere [33j . 

34j. In Fig. [T] we compare for (3 — 2 the 



We evaluated Eqs. (fT7f and (fT8f numerically using mathematica 
analytical results for characteristic functions and distributions with simulations obtained with an ensemble of 50000 
random matrices H of dimensions 200 x 200 from the GUE. The agreement is excellent. Unfortunately, there were no 
experimental data available for this case because complete violation of the T invariance could not be achieved. For 
8 = 1 we found that R s (k) is best evaluated using the Efetov variables 0q, 61,62 (0 < 6 < tt,0 < 6i t 2 < 00) [29|. 
These are related to the A's as Ao = cos do and A1.2 = cosh(6'i ± 62)- The numerical evaluation of the fourth derivative 
needed for the computation of the distributions P s (x s ) is not feasible. We therefore instead determined them with 
the help of Eq. ([7|l. considering a cut-off for k. This approach works well for a sufficiently flat distribution, whereas, 
if it is highly localized, it is advantageous to consider the corresponding characteristic function instead. 
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FIG. 1: Comparison of analytical and simulation results for /3 = 2: (a) Characteristic functions (b) Distributions of the real 
and the imaginary parts of S a b for the choice of parameters M — 3, E — 0.25, v = 1, 71 = 0.8, 72 = 1.0, 73 = 1.2, a = 1, b = 2. 



We tested our analytical results for /3 — 1 with experimental data. The latter stem from microwave measurements. 
To realize a chaotic scattering system, a microwave billiard with the shape of a tilted-stadium 2CH22J was chosen and 
the resonator modes were coupled to the exterior via two antennas attached to it. An ensemble of several chaotic 
systems was obtained by introducing a small scatterer into the microwave billiard and moving it to six different 
positions [36( . For the determination of the off-diagonal S'-matrix elements a vector network analyzer (VNA) coupled 
microwave power into the resonator via one antenna and determined magnitude and phase of the transmitted signal 
at the other one. The frequency range was chosen such that only the vertical component of the electric held strength 
was excited. Then the Helmholtz equation is mathematically equivalent to the Schrodinger equation of a quantum 
tilted stadium billiard. The dynamics of the corresponding classical billiard is chaotic 35j . The S'-matrix elements 



were measured in steps of 100 kHz in a range from 1-25 GHz. Disturbing effects due to the coaxial cables connecting 
the VNA to the microwave billiard were eliminated by calibration. Furthermore, the fluctuation properties of the 
S'-matrix elements were evaluated in frequency windows of 1 GHz in order to guarantee a negligible secular variation 
of the dissipation in the walls and of the coupling between the modes in the coaxial cables and the resonator modes, 
which are incorporated in the coupling vectors W c in Eq. @. More details concerning the experimental setup and 



the measurements are provided in [20j, |2l|. In Figs. [5] and [31 we test the analytical results with experimental data 



for the frequency ranges 10-11 GHz and 24-25 GHz, corresponding to a ratio of the average resonance width T and 
average resonance spacing d, T/d = 0.234 and, respectively, T/d = 1.21. The agreement is very good. 




FIG. 2: Characteristic functions Pi and P2 corresponding to the real and the imaginary parts of S12 for B = 1. Comparison 
between the analytical results and the microwave experiment data for the frequency range 10-11 GHz [20h22i ]. 




. (b) 




FIG. 3: Distributions of the real (Pi) and the imaginary (Pa) parts of Si 2 for B 
and microwave experiment data for the frequency range 24-25 GHz [2CH22I]. 



1. Comparison between analytical results 



To conclude, we solved the long-standing problem of deriving the full distribution of the off-diagonal S'-matrix 
elements for arbitrary values of the parameters involved. We accomplished this task by introducing a novel route to 
the sigma model based on the characteristic function. We verified our analytical results with numerical simulations 
and found excellent agreements. Furthermore, we tested our results with experimental data obtained with microwave 
billiards, and thus presented a new confirmation of the random matrix universality conjecture. 

This work was supported by the Deutsche Forschungsgemeinschaft within the Collaborative Research Centers 
SFB/TR12, "Symmetries and Universality in Mesoscopic Systems" and SFB634 (subproject "Quantum Chaos and 
Wave-Dynamical Chaos"). 



APPENDIX: EXPLICIT EXPRESSIONS FOR THE KiJ 



« u = -(9/8){p+m 1 / 2 } + , 

«2! = -(l/4)(128p°pg + 14p+P^ + 32p-p-) + {Ze^(p-q+- p -rt)}- + {e-*^-r 6 -} + , 

K22 = -(l/4){(p+p+ - 4g-r-)m} +! 

« 3 i = {2[(p+p+ + g-r-)m 1 ^ + 2 (8pK + p+p+ +p-p 6 ")/ 1/2 ] (e^" + e~ ia ^)}_ 
+{2[(p+^ + 4p-p+)mV 2 +p+pn 1 /2] (e"°V + e <a *r+)}_ + { [l6p°(2p° a p+ - 3p°p+) 

-6Pa (Qalb + r t r b) + 2 Pt( A( ii r i - 1a r a) ~ ^PaiP^Pb ~ 2 PaPt) + ^t (la % + r a r b ~ PtPb) 

+(e-***/2) g -(4p+r 6 - - 3p+«") + (e^/2)r-(4p+« 6 - - 3p+r-)]m x / 2 } + 

+ { [(e-^/2)q-{2e-^q-r- - 8e«*r+r+) + (e^ /2)r~ (2e^q^ r - Se^g^m 1 / 2 },, 

^32 - {p+ [{p + a P + a + 2q-r- a ) + (3/2)(e-^ff-«- + e^r")]™ 3 / 2 } , + {(2p+p+ + GY^e-**?- + e^r-)m 3 / 2 }_, 



Ml = 32[2p° Q p° a (p 6 - - e*+qt)te e'^r*) + 2pgpg(p~ + e' i2i '«+)(?£ + e«*r+) 
+ p°p°((p+ + e- i2 ^")(p+ - e-^rj") + (p+ + e i2 ^r")(p+ - e °V))] 
+ 266p2p2rfpg + (pt + e'^q-fipi e^r^f + ( P + a + e°*r") V - e°V) a 
+ 4[(p+ + e-^ 9 -)(p+ - e°V) - 2(p- + e-* 2 V)(p 6 - - e^q+)} 
x [(P+ + e^'r-)(p+ - e-^r-) - 2(p" + e**r+)fo - e"^r+)] , 

K42 = -32p°p° [(p+ + e-^q-M + e^r-)m + (p+ - e^V)^ - e _<2 V)J] 
-2[(p+ + e- i2 ^-)(p+ - e ia ^~) - 2(p" + e- ia ^+)(p fe - - e* 2 *^)] [(p+ + e^r") 2 ™ + (p+ - e^-) 2 /] 
-2 [(p+ + e^r-)(p+ - e"^r 6 -) - 2(p" + e^r+)(p- e"^r+)] [(p+ + e -°^-) 2 m + (p+ - e a ^)H] , 

-43 = (p+ + e-^q-) 2 (pi + e°*r-) a m a + (pf - e°V) a (rf - e^ 2 ^) 2 * 2 . 

In the above equations, we introduced the notation {£ (a, 6, Z, m, '0)}± := £ (a, 6, Z, rn., ip) ± £(6, a, m, Z, —ip) with £ an 
expression involving a, 6, Z, m, ip. 
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